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Abstract 

, We prove that the MaxweU-Klein-Gordon system on R^^* relative to 

^ ' the Coulomb gauge is locally well-posed for initial data in iJ"^"*"^ for all 

I £ > 0. This builds on previous work by Klainerman and Machedon Q 

^SJ . who proved the corresponding result, with the additional restriction of 

small-norm data, for a model problem obtained by ignoring the elliptic 
features of the system, as well as cubic terms. 

O ; 1 Introduction 

' The purpose of this paper is to prove local well-posedness (LWP) of the Maxwell- 

a Klein- Gordon (MKG) equations on ]R^+*, relative to the Coulomb gauge, for 
initial data in H^^^ , any e > 0. This result is optimal in the sense that the 
^ . critical Sobolev exponent for MKG on R^+^ is Sc = 1, and one does not expect 

' well-posedn ess in H'' for s below this critical value; see the introduction in 

. and section 1.3 below, where we also make some remarks on the open question 

' of well-posedness in the critical data norm H^. 

The analogous result for a hyperbolic model problem, obtained from the 
MKG system ^ below by setting the non-dynamical variable Aq = Q and 
ignoring all cubic terms, was proved by Klainerman-Machedon j^, for small- 
norm initial data. That result was reproved, using different norms, and without 
any smallness assumption on the data, in the recent survey article 0. The 
proof given there also used some ideas from ||^ , where the corresponding model 
problem for the Yang-Mills equation is considered. 

'Current address: Inst. f. Mathematik, Univ. Wien, Strudlhofgasse 4, A-1090 Wien, Austria 
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The present work builds on the treatment of the model problem in j?] : To 
obtain a priori estimates on solutions of MKG with the requisite regularity, 
we complement the bilinear estimates proved there with estimates for cubic 
terms, and terms involving the non-dynamical variable, which satisfies an elliptic 
equation. It should be emphasized that the difficulty is to obtain LWP when s 
is very close to Sc = 1- If s is sufficiently large, one can prove LWP by much 
simpler me thods than those employed here. See section 1.3 and Remark |^ in 



section 



1.4, 



Our method here can be modified^ to treat the full Yang-Mills system in 
proving LWP in but only for initial data with small norm. This 

extends the result of Klainerman-Tataru |^ on a model equation for Yang- 
Mills. The reason for the small-norm restriction is that the elliptic equation in 
the Yang-Mills system relative to the Coulomb gauge is far more complicated 
than the one for MKG, and not in general globally solvable. To avoid this 
problem one can include the elliptic variable in the Picard iteration. Then to 
close the iteration one must assume small-norm data, since there is no way of 
compensating for large data by letting the existence time go to zero, as one can 
for an iteration involving only hyperbolic equations in a subcritical regime. Of 
course, using Picard iteration for an elliptic equation seems somewhat contrived. 
A better approach for Yang-Mills on may be to work in the temporal gauge, 
as Tao |l^ has successfully done for the case of M^+'^. We hope to address this 
in a future paper. 

Most of the previous work on MKG has been in dimension 1 + 3. Let us 
summarize the known results for this case. LWP in the energy norm was 
proved by Klainerman and Machedon jij . By conservation of the MKG energy, 
their result implies global well-posedness. In particular, they recovered an earlier 
global regularity result of Eardley and Moncrief for smooth data. Cuccagna 
[0 proved LWP for small- norm data in i?*, s > 3/4. For 1-1-3 dimensions, the 
critical regularity is Sc = 1/2, but the question of LWP below s = 3/4 remains 
open. In both § and the Coulomb gauge is used. More recently, Tao |t7| 
has proved small- norm LWP for s > 3/4 using the temporal gauge, for the more 
general Yang-Mills equations. 

Our method here can be used to remove the small-norm restriction in the 
result of Cuccagna. The essential reason for this limitation in is that the 
elliptic variable was included in the iteration. If instead one solves the elliptic 
equation and reduces to a purely hyperbolic system, as we do here, this obstruc- 
tion is removed, and one can get a large data LWP result. A crucial fact needed 
to make this work is that in the Klainerman-Machedon bilinear estimates used 
by Cuccagna, the space-time derivative |£'t,a:| " acting on the product can be 
replaced by jD^;!"", as observed in [|ll[ (cf. also the remark in the Appendix), 
rendering unnecessary the decomposition in Fourier space used in j^]. See also 
Remark ^ in section |[ 

^We do not prove this here, but hope to address it in a separate paper dealing with the 
Yang-Mills equations on R-'^+'* in Coulomb as well as temporal gauge. Note that Yang-Mills 
essentially contains MKG as a special case. 
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Finally, we remark that our proof should generalize without difhculity to 
the higher dimensional case of MKG on with n > 5, giving LWP for 

s > Sc = In fact, the difficulty of the problem decreases with increasing 

dimension. 

1.1 The Maxwell-Klein-Gordon system 

The Klein-Gordon equation can be derived as a relativistic analogue of the 
Schrodinger equation for a free particle. It is obtained from the relativistic 
energy-momentum relation = p^c^ -I- m^c^, where E is the energy of the 
particle, m > its rest mass, p its momentum and c the light speed. Setting 
c = 1 from now on, and applying the quantum mechanical principle of replacing 
classical quantities by operators: 

Energy E — > i^, 

Momentum p — > iV, 

one obtains the free Klein-Gordon equation 

□0 = m^(j), (1) 

where (j){t, a;) G C and □ = df^d'^ — — + A is the wave operator on R^+" . Here 
we use relativistic coordinates t = x^, . . . , on the Minkowski spacetime 

with the metric diag(— 1, 1, . . . , 1); indices are raised and lowered relative 
to this metric, and the Einstein summation convention is in effect: roman indices 
j, k, . . . run from 1 to n, greek indices fi,!/, . . . from to n. We write for 

and dt — do. We shall use 5Rz and 3z to denote the real and imaginary 
parts of z £ C. 

The coupling of (|l|) to an electromagnetic field represented by a potential 
A^{t, a;) S R is achieved by the so-called minimal substitution 

9p — > D^^df, + lA^, 

where iA^^ acts as a multiplication operator. This gives 

Di^Df^cf) = m^(t). (2) 

which is the Klein-Gordon equation. It has an associated current density 

j^^^{mi:^)^^{0;^)-A^\.p\\ (3) 

satisfying the conservation law 

d^j^ - 0. (4) 

In fact, one has the general identity 9^3 [cf)Dt^(l)) — 3 (01)^13^0), so (^) follows 
immediately from (^. 
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The Maxwell-Klein-Gordon system is then obtained by coupling to the 
Maxwell equation 

d-F^,=j^ (5) 

where F^^, = df^Ai, — d^A^ is the electromagnetic field tensor and is the 
Klein-Gordon current density (||). The system (|5|),@,(||) is then what we — 
provisionally — call the Maxwell-Klein- Gordon system. We want to consider 
this as a system of second order PDE in the unknowns A^ and (/), but there 
is an obvious problem with this, since F^i, — and hence the observables, i.e., 
the electric and magnetic field vectors, whose components are entries of the 
matrix F^i, — are not uniquely determined by A^. This is known as the gauge 
ambiguity, and to resolve it one adds another equation to the system, a so-called 
gauge condition, which uniquely determines A^. The standard gauge conditions 
are (i) Lorentz: d'^Af^ = 0, (ii) Coulomb: d^Ai = and (iii) temporal: Aq = 0. 

In this paper, we shall rely on the Coulomb condition, which carries the 
advantage — as Klainerman and Machedon observed in [Q for the case of n = 3 
— that the bilinear terms involving derivatives turn out to be of null form type, 
and therefore have better regularity properties than generic products. Since 
the derivation of the null form structure in [|| uses the special vector calculus 
of n = 3, in particular the curl operator, we include a generalization of this 



argument to arbitrary dimension in section L5 



1.2 Main result 

If we add the Coulomb gauge condition Aj = to the MKG system (P), (^, 
(H) and expand, we get: 

AAo = -3(09^) + 101' Ao, (6a) 

□Aj = -'^{(fd;^ + 101' Aj - djdtAo, (6b) 

□0 = -2iA^dj(l) + 2iAodt(j) + i(5tAo)0 + Af'A^.ij) + m^cj), (6c) 

d^Aj ^ 0. (6d) 



In the rest of the paper, with the exception of section |1.3| , we will take n = 4. 
Thus, the unknowns are 



ol+4 



When convenient, we shall write A for the four- vector field (v4-')j=i_..._4. Initial 
data are specified at time t = 0: 

0L=o = ^^^ ^Mt^o = e (7b) 

where if** = {/ e ^'(R*) : (/ - A)''/^/ g L'^{K^)) and a, b are real vector fields. 
In view of the Coulomb condition ( |6d| ) , we must require 

d^aj = d^hj = 0. (8) 
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Observe that no data are specified for the non-dynamical variable Aq. This is 
quite natural, because Aq is determined by (j) and dt(f> at any time t by solving 
the elliptic equation (|6a|). 

Theorem 1. For all s > 1, the Cauchy problem on is locally 

well-posed. 

Local well-posedness here includes (a) existence of a local solution 

Aq eC{[0,T],H^) nC\[0,T],L^) (9a) 
Aj,(l)e C{[0,T],H'')nC^{[0,T],H'-^) (9b) 

up to a time T > depending continuously on the iJ^-norm of the initial data; 
(b) uniqueness of the solution; (c) continuous dependence on the data; and (d) 
persistence of higher regularity. A more precise statement, for an equivalent 
system, can be found in Theorem ^ section 1.5. In particular, the uniqueness is 



proved not in the class (M), but in a smaller space determined by the iteration 



norms; see (|20[ ). 

To prove Theorem |l| we shall in effect eliminate the nondynamical variable 
Aq from the equations, by solving the elliptic equations. This leaves us with a 
system of nonlinear wave equations, which we then prove is locally well-posed. 
Once this has been achieved, we can go back to the original system (|^), and 
conclude that this is also well-posed. 

Let us be more precise. We introduce a new variable Bq — OiAq. Applying 



to ( |6b| ) and using ( |6d|) yields 



ABq = {<l>d,(f,) + (101' Aj). (10) 

Now we eliminate Aq and dtAQ = Bq from (^) and (^) by solving (^^ and 
Thus Aq = Ao(0) and Bq — Bq{A, (p) are nonlinear operators. Since the 



Coulomb condition (6d) turns out to be automatically satisfied because of the 



constraint (jg), we obtain a system of nonlinear wave equations 

nA^M{A,(j)), (11a) 
□0 = A/'(A,0), (lib) 

where M. and are certain operators^, nonlocal in the space variable, which are 
sums of terms of the following types: (i) bilinear and higher order multilinear 
expressions involving A and (j) and their first derivatives, (ii) terms involving 
Ao((/>), and (iii) a linear term m^(/) in ( |llb| ). Moreover, all the bilinear terms have 
a null structure, due to the Coulomb gauge, and for these terms one already has 
good estimates (see @, and also ||] for the case of Yang-Mills; here we shall rely 
more particularly on variants of these estimates proved in 0). We complement 
these with estimates for the higher order multilinear terms and terms containing 
((/>), and local well-posedness of the system (^l]) then follows by the general 
theory developed in the author's paper [O. 



See section 1.5 for precise definitions 
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Then the original system (||) is also locally well-posed, by reversing the steps 
leading to (|ll|). That is, if {A, (j>) has the requisite regularity (see (pO|)) and solves 
( pll ) on a time-slab, and if we set Aq = Ao{cj)), then OiAq — Bo{A,(j)) in the 
sense of distributions and the triple {Aq, A, 4>) solves on the same time-slab. 

Thus, we show that the systems (||) and (^Tj) are equivalent for sufficiently 
regular solutions. 



1.3 Scaling, optimality and the null condition 

As for many other field theories, there are two types of "critical" behaviour 
associated to the MKG system on R^+". On the one hand, there is the critical 
regularity Sc such that the homogeneous initial data space H'^" is left invariant 
under the natural scaling transformation associated to MKG: 

(t , x) , 0(t , x) — > XAf, (At, Aa;) , A0(Ai , Aa;) , (12) 

where A is a positive parameter.^ Since 

||A/(A.x)||^.=A^-("-2)/2 11/11^., (13) 

we conclude that Sc = ^-^^ ■ In general^ one expects field theories to be locally 
well-posed (LWP) for s > Sc and ill-posed for s < Sc, we say more about this 
below. In the critical case s = Sc one expects some type of weakened well- 
posednessj^ for data with small norm. 

On the other hand, there is the energy-critical dimension n such that the 
critical regularity is at the level of the energy:^ Sc — I- For MKG this means 
n = 4, which is the dimension we consider in this paper. For field theories in 
general, one expects global regularity in the critical dimension, as well as in 
subcritical dimensions (sc < 1), and breakdown of regularity for large data in 
supercritical dimensions (sc > 1). 

As mentioned above, the global regularity is known in the subcritical di- 
mension n — 3 for MKG, but the question of global regularity in the critical 
dimension n = 4, even for data with small energy, remains open. By conser- 
vation of energy, a LWP result, for small- norm data, at the critical regularity 
Sc = 1 would settle this question in the affirmative, but it is perhaps more 
realistic to expect a more direct proof of global regularity in analogy with the 
results of Tao |jl^, |l^ for wave maps into a sphere. It is to be hoped that our 
almost optimal LWP result will play some role in any such result. 

The expectation of ill-posedness for s < is based on the scaling ( p^ ) and 
13[). First, if blow-up occurs for smooth, compactly supported data, then one 



•^By this we mean that if A^,(/) solve MKG, then so do the rescaled fields, although the 
rest mass changes from m to Am. 

*See wA Section 1.3] for further discussion and references. 

^For example, one does not expect smooth dependence on initial da ta, which rules out 
proof by iteration. A good example is wave maps into a sphere; see Tao fllSl Section 1] for a 
summary of the regularity results for wave maps. 

^MKG has a conserved energy which is at the level of the data norm; see M]. 
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can construct data in H^, s < Sc, with arbitrarily small norm, for which there 
is no local existence; see, e.g., pp 98-99] for this argument. However, this 
is not a very convincing point to make here, as we do expect global regularity 
for MKG on R-'^+'*. We can show, however, that it is impossible to prove any 
well-posedness result for s < Sc using an iteration argument based on estimates. 
The idea can be illustrated by the following example: As is well-known, the 
algebra inequality 



(14) 



holds for H'^ (R") iff s > n/2. A rather crude way of ruling out the range s < n/2 
is to observe that if (|lj) holds, then by rescalin^ x ^ Xx and letting A — > cxo, 
we get 1 < A'-"/2. 

This idea is easily applied to the iteration for MKG written in the form ( [ll|) . 
Let us take m — here to make the system scale invariant. If a = 6 = in (iTaj) 



and 01 = in (jTy), then the first iterate of A solves| DA^^) = -PSy(0(o) V0(o)) 
with zero data, where 0^°) is the solution of D^'"' = with data (0o,O). If we 
can prove LWP in by iteration, there must be an estimate 



sup 

0<t<l 



< 



1100 



(15) 



for all ipo with sufficiently small norm. Now assume s < Sc- We then claim that 
( p^ implies A^^^ = 0, which is absurd. Indeed, given T > 0, apply ( p^ ) to the 
rescaled iterate _ 

A'^^Ht,x) = XA^^\Xt,Xx) 

at time t = T/X. As A ^ oo this gives 



A^ 



< 



(A^ 



whence A'^^^T) = 0. 

Remark. This argument has nothing to do with the null condition, of course. A 
more careful analysis (see [Q, Section 1] ) suggests that for a generic equation of 
the form 

□u — udu 

on one needs s > max(2Y^, '^^^) in order for the iterates to stay in H^, 

and this is consistent with Lindblad's counterexamples llffl. However, if the 
right hand side is replaced by a null form expression like ( p3| ) or ( p^ , then one 
only needs s > max(^2^ ilzlL)^ go the null condition improves matters when 
n < 4. 



As remarked already, the main difficulty is to prove LWP when s is very 
close to Sc, whereas simpler arguments can be used for larger s. Let us be more 

^In the limit A oo, the inhomogeneous Sobolev norm scales like iif*. 



'Here V denotes the projection onto divergence free vector fields. See section 1.5 
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precise. Observe that relative to Lorentz gauge, MKG on K-'^+" is a system of 
nonlinear wave equations of the schematic form (see |Q, Section 1] ) 

Du^udu + u^, (16) 

and for this system LWP for s > n/2 can be proved by standard methods, 
just using the energy inequality for the wave equation and Sobolev embeddings. 
This can easily be improved to s > ^^^^ by using a L'^L'^ spacetime estimate 
instead of just Sobolev embedding. For n = 4 this gives LWP for s > 3/2, 
which is still one quarter of a derivative above what one expects (cf. remark 
above) from the analysis of the first iterate of ([T6|), namely s > 5/4. No proof 
of LWP of (|l6|) in this range seems to exist in the literature, but it should be 
obtainable using the spaces iJ*'^ (see section and bilinear estimates for 
the homogeneous wave equation of the type first proved in However, to go 
below the regularity 5/4, one really needs the null condition, which seems to 
rule out Lorentz gauge. Of course, once a LWP result has been proved in one 
gauge, one can in principle use gauge transformations (see Q) to transfer this 
result to other gauges; but to make this rigorous requires sufficient regularity of 
the solutions, and we will not consider this question here. 

1.4 Function spaces 

Here we define the spaces that we make use of. See for more details. 

The Fourier transform of f{x) [resp. u{t,x)] is denoted /(^) = J-f{C) [resp. 

We say that a norm ||-|j, on some space X of tempered distributions, depends 
only on the size of the Fourier transform if 

l^^l < \v\ =^ \\u\\ < \\v\\. 

(Here we assume, of course, that the Fourier transform of any element of X is 
a function.) 

If X and y are two normed spaces, the notation X '-^ y means continuous 
inclusion. 

For any a e M we define Fourier multiplier operators A", A^j: and A" by 

A^/(0 = (i + iei')"^7(o, 

A^u{r,0^il + r' + \^\Y'u{r,0, 

A^m(t, 0=1 + ^ Hr, 0- 

\ l + r2 + |e|7 

It should be remarked that the weight of A" is comparable to(l+||r| — |^||)", 
but the former has the advantage of being smooth. 

The Sobolev and "Wave Sobolev" spaces and iJ*'^ are given by the 
weighted norms 

II/IIh^- = I|A'/IIl2(r4) and = \\A'Atu\\^,^^,^,^ . 
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We shall also use the related space H.'^'^ defined by 

WAh-.o = WAh-.o + ||<9tu||^.-i,e - ||A"~^A+A%||^2 . 

In view of Plancherel's theorem, these norms depend only on the size of the 
Fourier transform. It is an important fact that when > 1/2, the spaces H'^'^ 
and 7i*'^ can be localized in time, since then the cmbeddings 

H'^'^ ^Cb{R,H'') and W''' ^ Cb{R, H") D Cl{R, H'-^) (17) 

hold. See g, Section 3]. 

Since i2(|^|2 c Ll^iR-^) C S'{R^), we may define 

H^=T-'[L\\^\'dO]. 

Thus is a Hilbert space with norm = J^4 |^|^ |/(C)|^'^C- We remark 

that if ~ {/ : V/ G L^}, then is obtained by identifying elements of 
W-^ differing by a constant. Observe also that S is dense in L'^{\^\ d£^), hence 
in H^. We shall use frequently the fact that 

^ L\R'^). (18) 

In other words, ||/||i4_< ll/ll/fi- This holds by the Hardy-Littlewood-Sobolev 
inequality (see Stein |lj. Chapter V]). 

If <Y is a separable Banach space of functions on R**, and 1 < p < oo, we 
denote by L^{X) the space L^^R^X) of -^-valued functions. In particular, we 
write 



\W\\mL%) = (^^II"(*>-)IIl,(r4) 



with the usual modification if p = oo. 

We also need a version of this last norm which only depends on the size of 
the Fourier transform: If m £ iS' and w is a tempered function, set 

WAc'^id) = '^'"^P {/i+i \u{T,£.)\^{'r,CldTdi ■.veS,v>0, WvW^p'^^^'^ = l| , 

where ^ — p ~^ y ^^"^ ^ ^ q ^ V' ^ti^x) the corresponding subspace 
of S'. Then |H|£P(£9) is a translation invariant norm on £((£^). Note that 

CUCl)^L^{R'+h^^d 

\\u\\c''{c'') — \\Al^{l'') whenever u>0. (19) 

We refer the reader to Section 4] for more details on these spaces. 

We can now make precise the regularity statement (^. The solutions we 
obtain arc in the following spaces: 

Aa e C{[0,T],H^) nC^{[0,T],L^), (20a) 

A.eH^^'nA-^AZHcliCD], (20b) 
e (20c) 
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where > ^ and 7 > depend on s. For technical reasons, it is useful to iterate 
Aj and 4> in these global spaces, but in the end we are only interested in their 
values on a time interval [0, T] whose size depends on the norms of the data. 
Since the space 7i^'^ can be localized in time, this presents no problems. 



Remark 1. The auxiliary space C\{C^) in (20b) is necessary when s < 5/4. See 
Theorem 8.2 in 0| and the remark following it. 

Note: Throughout the paper, we use the convenient shorthand < for < up to 
a positive multiplicative constant C. Usually C is completely innocuous, and 
only depends on parameters that may be considered fixed. There are exceptions, 
notably for Lipschitz estimates (then C is only "locally" constant), but these 
are clearly pointed out. 



1.5 Reformulation of the MKG system 



As discussed in section 1.2, an important step in our proof is to recast the MKG 
system (^ as a system of nonlinear wave equations (^l|). Here we describe this 
in detail. 

As was shown in the first terms on the right hand sides of equations (6b) 
and (|6c|) can be expressed, due to the Coulomb condition (6d), in terms of the 
bilinear null forms 



Qjkiu,v) = djudkV - dkudjV. (21) 

Since the argument in Q was special to the case n = 3, we include here a 
proof of this fact which works for any dimension. First, let V be the projection 
onto the divergence free vector fields on M"*. In terms of the Riesz transforms 

R, = (-A)-^a„ 

VXj ~ Xj + RjR^Xk = R''{RjXk — RkXj). 

Observe that V is bounded on every L^, 1 < p < 00, since this is true for the 
Riesz transforms (see Stein [Q). Moreover, it is clear that the Riesz transforms, 
and hence "P, are bounded on any space whose norm only depends on the size 
of the Fourier transform, in particular on any Sobolev space . 

Since dj{udkv) — dk{udjv) = Qjk{u,v), it follows immediately from the def- 
inition of V that 

Viudjv) = R''i-A)-iQjkiu,v), (22) 

whence 

V{-^0~^) = 2i?'=(-A)-5Q^fe(sR0,30). (23) 

Also, 

2djuVX^ = (-A)-^ [WX^ - R^X^]) , 
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as one can see by expanding the right hand side. Therefore, if A is divergence 
free, so that VA = A, then 

2A^dj(b = Qjk (0, (-A)-^ [R^A'' - R'^A^ ) . (24) 



Remark 2. The calculations leading to the identity ( |22| ) are certainly justified 
when u and v belong to the Schwartz class iS(R^). Moreover, both sides of the 
identity are bounded bilinear operators of {u,v) G H" x into H~^, where 
s > 1. Thus the identity holds for all u,v £ -ff*, and we conclude that ( |2^ ) holds 
for all (j) with the regularity (pO^ ), since by (|l^) this implies (j) £ Cb{R, -ff"). To 
bound the left hand side of (|22|), use first the dual 

ll(-^)"^/|L=(E^) ^ II/IIl^/3(e^) (25) 
of (|l8|). Since V is bounded on L^, it then suffices to observe that 

ll'«<9jw||^4/3 < ||m|Il4 \\dju\\L^ < hWm \Hm > (26) 



where we used (|18|). To prove boundedness of the right hand side of (g2|), it is 
enough to show 

||(-A)-^(/g)||^_, <|l/||^._,|l.g||^._,. 

This can be reduced, via the self-duality of L^, Plancherel's theorem, and the 
Cauchy-Schwarz inequality, to the fact that |^|~^ (1 + belongs to 

L^(M^), since s > 1. Similar, but simpler, considerations show that the re- 
maining bilinear and cubic terms in (||a,b,c) and (]28|a,c,d) are bounded into 
C6(R, L*/'^(R*)) when regarded as operators on Aq, A, (j) in the class (pO|). For 
example, for a cubic expression uvw we have by Holder's inequality and ( p^ 
that 



uvw 



Il4/3 < 111^1^4 \\V\\^, \\w\\^, < \\u\\^^ \\v\\^^ \\w\\J,^ . (27) 



Returning to the main thread of our argument, we now use the null form 
identities derived above to arrive at an equivalent formulation of MKG: 

AAo - -3((^5^) + Ao, (28a) 

AdtAo - -W^ {(fd~4>) +d^{\(t)\^ Aj) (28b) 

UAj = 2i?^-(-A)-^Q,fe(K0,30) +P(|(/)|' A,) (28c) 

□</, = -iQ^-fe (^0, (-A)-^ [WA'' - R^A^]^ (28d) 
+ 2iAodt(j) + i{dtAo)<j) + A^'Af,<j) + rn^cj). 
This system acts as a stepping stone between (^ and (pT|). 

Proposition 1. The systems (|^) and ( |28|) are equivalent. More precisely, any 
local solution of (|^) with the regularity (|20|) and divergence free initial data is 
a solution of (E81) and vice versa. 
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Proof. To go from to ([2^), observe that Aj is divergence free by (6d); apply 
to ( ph|) t o get equation (|28b ); apply V to (6b) and use ( |23| ) to get ( |28cD ; 

finally, (p8d|) follows from (|^^sing (|2|). 

To go the other way, observe that by (|2^), the right hand side of ( p8c| ) is 

divergence free; thus Od^Aj = 0, and since the initial data of Aj are divergence 



free, (6d) follows. Then, in view of (p4|), (p3q) and (28d) are equivalent. Finally, 



to go from (28c) to (6b), it suffices to chec k th at the right hand side of the latter 

□ 



is divergence free. But this follows from ( p8q ). 

Once the system has been written in the form (|2^) it is easy to eliminate Aq 
and dtAo and obtain the system of wave equations (|l^). We now describe this 
in more detail. 



Lemma 1. Given <j) in the class (|9q), equation (28a) has a unique solution 
Aq £ on every time-slice {t} x and these solutions assemble to a space- 
time function Aq = Ao{4>) e Cb(R, H^). Moreover, we have bounds, on every 
time-slice {t} x R"*, 

\\AoUi<^\\Mh^ 

and 

||Ao(0) - Ao(V)|Ihi < U - + - Mh^ , 

where the suppressed constant depends polynomially on ||(/>l|^i , l|'01l_f/i j llf^i'/'llis 
and \\dtip\\L^j but is independent oft. 

This is proved in section ^. 

Next we consider (28b), with dtAo replaced by the new variable Bq: 



ABo = -W{(l}dj(f,)+d^{\(f>\^Aj). 



(29) 



Lemma 2. Given {A,(jji) in the class (|9b|), the equation ( p9| ) has a unique so- 
lution Bq E on every time-slice {t} x , given by 



Bo - i?^(-A)-5 [3(09j0) - 101' A, 



(30) 



and the solutions assemble to a space-time function Bq — B(){A, (j)) G Ch 
Moreover, we have bounds, on every time-slice {t} x M*, 

\\B4l^ <C{l + \\A\\H^)m\^ 

for a constant C independent of t, and 

\\Bo{A,cj>) ^ Bq{A' ,^')\\^. <\\A- A'W^, + U - cf^'W^, , 

where the suppressed constant depends polynomially on ||^||^i , ||v4'||^i , Ht/'ll^i 
and \\4>'\\h'^, but is independent oft. 



Proof. To see that ( |30| ) is in L^, first apply (|25|), then estimate as in ( Pq ) and 
(|7|). That (|0|) is the only solution can be seen by taking the Fourier 
transform of both sides of ([29|). □ 
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In view of the above lemmas, (^8|) implies (jllj), with 

M = {Mi,...,Mi), Mj ^M,,i+Mj^2, AA^AAi + .-.+A/'e, 

where 

M2 - 2^0(0)5*0, 
M^iBa(A, 4,)<j,, 

Mi = -[Ao(0)]'0, 
M, = \A\^ct>, 
Mq = rri^ij), 



and 1^1^ — AjA' in the next to last line. 

Arguing as in Remark H, and using Lemmas |^ and ^, it is readily checked that 
the multilinear expressions in and M are all continuous maps into Cfc(]R, L'^/'^) 
[or Cb(R, i/^^) in the case of A^j.i] for (v4, in the class (pO|b,c). However, 
proving the following theorem requires much more sophisticated estimates. 



Theorem 2. The system of wave equations (|ll[), with A4 and M defined as 
above, is locally well-posed for initial data in H'^ , all s > 1, in the following 
sense (all pairs {A,(j>) are understood to belong to the class ([20^3, c) in what 
follows): 

(a) (Local existence) For all initial data (Q) there exists a T > 0, which 
depends continuously on the norms of the data, and there exists a pair 
(A, (f)) which solves ( p^ ) in the sense of distributions on (0, T) x R"* and 
satisfies the given initial condition. 

(b) (Uniqueness) If T > and we have two solutions {A,(j)) and {A',(j)') of 
(11) on (0,T) X with identical initial data, then they agree on the entire 
time-slab. 

(c) (Continuous dependence on initial data) If, for some T > 0, (A, 
solves (11) on (0,T) x M.* with initial data (|^), then for all initial data 
{a' ,b' , (f>Q, such that 

S^\\a- a'll^. + ||5 - + ||0o - ^'oWhs + Ui - Mh^~^ 

is sufficiently small, there is a solution {A',(j)') on the same time-slab and 
with these initial data. Moreover, we have 

\\A - A'W^, + WdtA - dtA'\\^^_, + 110 - + Wdtcl) - < CS 

uniformly in < t < T. 
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(d) (Persistence of higher regularity) Ifk is a positive integer and {A,(f>) 
solves ( pH ) on (0, T) x M"* with initial data in H"^'' (that is, holds with 
s replaced by s + k), then 

A,(l)e C{[0,T], H'+'') n {[0,T], H'+''-^) . 

(e) (Classical solutions) If the data belong to H'^^*' for every k, then the 
solution is smooth: 



A,4>eC°°{[{),T] 



X 



The proof of this theorem wiU occupy us in the next two sections. 

Here we want to show that Theorem ^ can be deduced from Theorem ^ It 
clearly suffices to demonstrate the equivalence of the systems (||) and (|Tl|). The 
remainder of this section is devoted to a proof of this fact, assuming that the 
conclusions of Theorem |^ hold. 

Proposition 2. The systems (|^) and ( pT| ) are equivalent for local solutions in 
the regularity class (|20|), with divergence free initial data. 

In view of Proposition [|, it suffices to show the equivalence of ( p8| ) and 
l|). We have seen already that (28) implies ([Til). The converse is not quite so 



obvious, but for sufficiently regular solutions it follows by some straightforward 
calulations and the fact, proved in section ^, that the only solution of the 
elliptic equation Au — |(/)|^ u is u = 0. For general H'^ data we then choose an 
approximating sequence of sufficiently regular data, use the persistence of higher 
regularity and continuous dependence on initial data, which hold by virtue of 
Theorem ||, and pass to the limit. 
We now turn to the details. 

Assume that (A, (j)) is in the class (pO[b,c) and solves ( pT| ) on a time-slab 
St = (0,T) X R"*, with initial data satisfying (0) and (|). Set = ^o('^)- 



Then (|28| ) is satisfied, but with dtAg replaced by Bq — Bo{A,4>) in ( 28b ) and 
( ^8d| ). Thus, all we have to prove is that the distributional derivative dtAo agrees 
with Bq on St- At first glance one may think that this is simply a matter of 
taking a time derivative of ( ^8a| ) and using the conservation law (jj) to conclude 
that Ai^t^lo = ^Bq, but this is a circular argument since the derivation of (||) 
is not valid unless we know that dtAg = Bq. 

In what follows, keep in mind that and Bq are real-valued. Applying dt 



to ( p8aD gives 

AdtAo = -5(0^0) + 25?(09^) Ao + |0|' dtAo. (31) 



Since (28c) and (H) hold, it follows as in the proof of Proposition!^ ^^^^ ^ 



divergence free. Therefore, (|2J) holds, and since (28d) holds (with dtAg replaced 
by Bo), we conclude that 

-d^(f) + A(/« = 00 = -2iA^dj4> + 2iAodt(f> + iBo(f> + A'^Af.cj) + m'^cf). 



14 



Using this expression for df(t) gives, after some calculation, 
Since 

-Qd^{4'd~^) =ABo-d^{\^fAj), 

we get 

-9(0^) ^ ABo - {\(f>\'^ Aj) + 2^(f>d~^A^ - 2^(fd^)Ao - Bo |</)|^ , 



Inserting this in ( |3l[ ) gives 

AdtAo = ABo-dm(f>\^A,)+2^{(f>d;^)A^-Bo\(l)\^ + \<l)\^dtAo. 

But 

A^.) ^ 2^(j)dj(t))A^ + \(t)\^d^Aj = 2^{(l)dj(l))A^ 

since A is divergence free, and so we finally get 

A{dtAo-Bo) = \cl)\^{dtAo~Bo). 
The above manipulations are justified provided 

dtAoeC{[0,T],H^). (32) 

If, moreover, 

Bo eC([0,T],iji), ^33) 

then it follows by the uniqueness result alluded to above (see Lemma ^ in section 
|) that dtAo = Bo in [0,T] x R"*, 

But (^2|) and (^3|) certainly hold under the additional assumption that the 
initial data (0) of A and (j) belong to H'+'' for every positive integer k. Leaving 
aside the proof of this assertion for the moment, we note that any / e H'' can 
be approximated in the i?** norm by a sequence belonging to every H'^^'', by 
convolution with a approximation of the identity, and if / is divergence 
free, then so is the approximating sequence. Combining these facts with the 
continuous dependence of A and (f> on their initial data (Theorem ||), and the 
continuity of the operators Aq and Bq (Lemmas and ||) , we conclude by passing 
to the limit that the equality dtAo — Bq holds in the sense of distributions on 
(0,T) X R'^ for all initial data^ satisfying (|). 

It remains to prove that 1^2) and (^3|) hold if the initial data (^ of A and (j) 
belong to for every positive integer k. For Aq, this follows by persistence 

of higher regularity (part (d) of Theorem ||), the inductive regularity step ( [73| ) 
in section 3.2 and Lemma ||in the same section. As for Bo, in view of ( |30| ) it is 
clear that, on every time-slice, 

and by Holder's inequality and Sobolev embedding it is easy to see that the 
right hand side is dominated by ||0||/f2 + ||^||^2 ||A|j^2. But if A and 4> 

have initial data in H^^^, then by persistence of higher regularity (part (d) of 
Theorem ||) wc know that A, </) G C([0, T],H'^). 
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2 Proof of Theorem 



Here we discuss the estimates needed to prove local well-posedness of the system 
Tlj), with M and Af defined as in section |1.5| . 

The local existence for the system (O) is proved by Picard iteration in the 



spaces (|20b|) and (20c), which are defined using the spacetime Fourier transform, 
and hence are global. However, since they embed in (|9b|), they can easily be 
localized in time. In fact, this time localization smooths out the singularity 
of the inverse of the wave operator, and — if done with sufficient care 
— allows one to handle large initial data by taking a sufficiently small time 
interval. These matters are considered in detail in the author's paper and 
also in Section 5], and we refer the interested reader there. 

Fix 1 < s < 2. (For larger s, the result can be proved by simpler arguments.) 
Let > ^ and 7, e > 0; these quantities depend on the choice of s, and will be 
specified later. Now define 



with norms 



All these spaces are complete (see Proposition 4.2]), and by ||^, Proposition 
5.6], Xi and X2 satisfy the hypotheses of |12, Theorem 1]. Consequently, by 
Theorem 2], the system ( pT| ) is locally well-posed for iJ* data if the following 
Lipschitz conditions^ hold: 

\\M{A, <t>) - M{A', cl^')\\y^ < \\A A'y^ + U - 4>'y, , (34a) 
\mA,^)-MiA',cj^')\\y^ < \\A-A'\\^^ + ||0-<^'1U,, (34b) 

where the suppressed constants depend continuously on 

IIA'IU^, 11011^^ and . 

In fact, these estimates guarantee that the conclusions (a,b,c) of Theorem || 
hold. In the next section we show how to prove parts (d) and (e) of the same 
theorem. 

It suffices to prove (^) with A4 replaced by A4j,k and with J\f replaced 
by A/i, . . . , A/5. Furthermore, in view of the multilinear structure, it suffices to 



^Keep in mind that M and J\f vanish at the origin, so if we take A' = and <j)' = 0, we 
simply get bounds for A4{A, <j)) and Af{A, 0). 
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prove (concerning the suppressed constants, see note below): 



||A^,,ilL, < 

IIA/-2 
IIA/-4 

IIA/'e 

||Ao(0)-Ao(<^')IIl,c.(^,)< 



y2 
y2 



< 



< 



< 



< 



< 



IU2 ' 



A\U II 
So(A0)IU, ll<^IU,, 



(35) 
(36) 
(37) 
(38) 
(39) 
(40) 
(41) 
(42) 
(43) 
(44) 



where Zi and Z2 are certain intermediate spaces, to be specified later, such that 

ll^ob,<||<^||^, + 11011^,, (45) 
||Ao(^)-Ao(0')b, <ll'^-'/''IU,, (46) 
\\B4z,<{l + \\A\\;,Jm%, (47) 
\\Bo{A, 0) - So(A', 0')IU, < P - ^'lU, + II'/' - 0'll;,, • (48) 

It should be emphasized that in the Lipschitz estimates (^J), (|4^ ) and (^8|), the 
suppressed constant depends polynomially on the norms and ll^'ll;^^' 

in the case of ( ^8|) also on and ||^'||^^. Observe that the estimate ( ^2| ) 

for the linear term is trivial, since the norms only depend on the size of the 
Fourier transform. 

The following was proved in ||^, Theorem 8.6]. 

Theorem. The estimates (Bq) and (|37|) hold provided 



1 



- < « < mm 1 -,- 



3 s 



< e < - min ( - 

4 

1 

2 ' 



7 



3£. 



(49a) 
(49b) 
(49c) 



Having fixed 9 and e satisfying these requirements, we define p and r by 

p 2 r 
and we choose q so large that 



1 - 61 - 2e, 



(50) 



4 

— < mm 

q 



29 



1 

1,1- - 
P 



(51) 
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Observe that as s — ^ 1, the triple {p, q, r) — > (1, oo, 2). Now set 

Poll, =||A^-Mn" 



IS, 



oWz. 



(52) 
(53) 



For easy reference, we hst here some estimates that we shall use (here p, q, r 
are defined as above): 



\A'-\~A)-\uv)\ 





< 


WWhs.b \\v 










< 












< 










ll"llL?''(Lf) 


< 






- 2e - 2s < 


-<2e, 

/? - 


u||^o,e+e-i 


< 












< 










11/511^^. 


< 


lA^/lliP, 


|5|Il,i 4 


' l|A^5lliP. 


l/IL,. , 



(54) 
(55) 
(56) 

(57) 
(58) 
(59) 
(60) 



where in the last inequality, 
fj > 0, 



1 

Pk 



1 



2 <pk < oo. 



we 



The inequality (^4|) follows from a theorem of Klainerman-Tataru 
give the details in an appendix. 

The Strichartz type estimates (|5^-p7[) are special cases of Theorem D]. 
(The [non-optimal] upper bound for 8/(3 in ( p7| ) guarantees that the pair {2p, (3) 
is wave admissible; the lower bound is chosen so that we do not exceed s space 
derivatives on the right hand side.) 

The inequality (58) can either be proved directly, using Plancherel's theorem. 
Holder's inequality, Minkowski's integral inequality and the Hausdorff- Young 
inequality, or it can be proved by interpolation, as in [Q, Section 6(vii)]). 

Inequality (^ is a special case of |0, Proposition 4.8]. 

The calculus inequality (|60| ) is Lemma 1 in Ponce-Sideris |^. 

As mentioned already, ([35|) an d (^7|) h old by @, Theorem 8.6]. We now prove 
the remaining estimates (t?6|) and (|38H48|) , thereby concluding the proof of parts 
(a,b,c) of Theorem ||. 

2.1 Proof of (H) 

Since the norm only depends on the size of the Fourier transform, we can ignore 
the projection V. More accurately, 



< 



'A 



JWyi 
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Thus, it suffices to prove 

\\A+^Al-\uvw)\\_^^ < \\u\\^.,e IH^s.e \\w\\^,,e , 

or, equivalently, 

]\uvw\\ij,^i,g+,-i < \\u\\^,,e \\v\\^s,g \\w\\^s,e , 
II -1 e-- II 

||A''A+ A_ '(wwti')||£i(£8) ^ hWn^.e IH-H^.e \\w\\^,,e . 

Since all the norms depend only on the size of the Fourier transform, we may 
assume that m, u, w have non-negative Fourier transforms, and we see that it is 
sufficient to prove (note that 7 + 2e < s — 1 by (^)) 

\\uvw\\fjo,e+,-i < ||M||^i,f) \\v\\^s,a \\w\\^,,e , (61) 
\\A-'^AZ'^^'{uvw)\\^,^^,^ < \\v\\^s,e \\w\\^,,e ■ (62) 

By (|5^) and Holder's inequality, 

\\UVW\\ ^0,0^,-1 < ||M||ioo(i4) ||w|Il2p(^8) ||w|ll,2p(^8) , 



and (|6l|) follows by Sobolev embedding and (57). 
Using (|9|) and (|l|), we get 

\\A-^AZ"'^{UVW)\\^,^^,^ < \\uvw\\^ii^^2) < I|w|Il~(L4) lkllL?(L|) M\Lj{Ll) 

Now use Sobolev embedding and (|55|). 

2.2 Proof of (H) 

We have to show 

By M) and {M), 



||uu||^,_i,e+,_i < IjA" ^(uu)||^p^^2^ 



< A 



The desired estimate now follows by Sobolev embedding, since | < s — 1. 

2.3 Proof of (H) 

We must prove 
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By H) and ®, 



Now apply ( |55| ) and (pq). Note also that ||u|j^^^^8/3^ < [[A'* -^u||^r^^8/3p since 

A^"^ is bounded on for all 1 < p < oo and 5 > 0. In fact, A^"^ corresponds 
to convolution with an function; see Stein p| . 

2.4 Proof of (13) 

It suffices to show 

By H) and ®, 



Now apply Sobolev embedding, and use (Fll). 



< llA'*-i- 



2.5 Proof of (0) 

It suffices to show 

\\uvw\\Hs-i.e+,-i < \\w\\^,,e 

but this was proved above; see the proof of (pq). 



2.6 Proof of (P) and (g4D 

These follow from Lemma |^, which is proved in section |^. 



2.7 Proof of (H) and (H) 

Since 



it suffices, taking into account the multilinearity of the terms inside the brackets, 
as well as the estimates (^) and (^4|), to show that 
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The former is exactly (54), and the left hand side of the latter is < 

\\{-A)-^{uvw)\\^^i^^,^ + \\{-A)^{uvw)\\^p^^^,y (64) 

Here we applied the following useful result, which is an immediate consequence 
of Lemma 2(ii) in Chapter V of Stein fl^ . 

Lemma 3. For a > and I < p < oo, 

iiA"/iL.<ii/iL. + ll(-Ar/viL„ 

where the suppressed constant only depends on a. 

Returning to the sum (|6^), note that by Sobolev embedding, it is < 



where 



1111 „/l 1 
— = - + - = - + 2 - + — ,, 
ai 2 g 4 V8 2qJ' 

1 3-s 1 1 /2-s 1 
2 



a2 4 q 4 \ 8 2q 



Thus 



ll^^llL^i:'') ^ II^IIl^^cl^) II'^IIl^'cl^) IMLriLt) 

where 

1 1 1 1 2-s 1 



1,2 



(3i 8 2g' P2 8 2q 
Using ( [49b| ) and (|l]) it is easily checked that 

5 8 8 

- + 9 + 2e-2s<2-s< — < — <2e, 

2 fh Pi 

so we may apply (^^ to finish the proof. 
2.8 Proof of (0) and (j^ 

We prove (^7|); the same proof gives ( ^8|) if one exploits the multilinearity of the 
terms defining Bq. 

First observe that by Lemma |3[ 

W^^'^MlmLl'-') - 11^011^.(^3/3) + ||(-A) — Bo||^^(^8/3). 

Therefore, by Sobolev embedding, we have to estimate 
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where 

1 5 1 5 s - 1 



Qfl 8 Q!2 8 4 

Since Bq is given by (|30|), and since the Riesz transforms Rj arc bounded on 
L^, 1 < p < oo, we see that it is enough to prove 

ll™^i'llL;-(ix'=) - IMh"." Ikllff-.f Ikllff^.f ■ 
By Holder's inequahty, 

IIw^^w'IIl-cl"^-) < ' 



where 



1 ls-1 1 



/32 2 4 ' 72 4 4 
Now apply ( p6| ) and Sobolev embedding. 

3 Higher regularity 

Here we prove parts (d) and (e) of Theorem 
3.1 The persistence property 

The key to proving part (d) of Theorem || is to establish, for fc = 0, 1, 2, . . . , 

\\A''M{A,^)\\y^ <afe{||A'=A||^^ + ||AV|Uj+/3fc, (65a) 
||A*'-A/-(A,0)||^^ <afe{||A'=A||_^^ + ||AV||;,J+/3fc, (65b) 

where 

• Q!fc depends continuously on and ll'^ll;^'^' 

• /3o = 0, 

• /3fe, for fc > 1, depends continuously on || A'^^-'^A||^ and ||A'^^-'^(/)||^ . 



The case fc = is of course true by (|3J), but it is useful to include it here for 
technical reasons. 

In the absence of the lower order term Pk , we could now appeal directly to 
[HI Theorem 2], to conclude that part (d) of Theorem || holds. However, we 
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can easily modify the proof given in to cover this more general case, as we 
demonstrate below. 

First, however, let us dispose of proof of the above estimates. Observe that 
we have the equivalence of norms 

|a|<fc 

This is trivial in view of the fact that the norms only depend on the size of 
the Fourier transform. It is therefore clear, from the multilinear structure of 
^A and Af, and the product rule for derivatives, that (^) follows from the very 
estimates proved in section ^ The only exception is the nonlinear operator 
Ao{(p), for which we need the following estimate, replacing (^3|): 

Lemma 4. If A'^cj) e X2, then 

\\d:M\Lrim)<-/kmxJ\\^'4x.+M\\^'"'4xJ f^'^'^l^ l"l 
where 7^ and rjk are continuous functions. 



This is proved in section 4.3 



Let us now turn to the proof of Theorem g, part (d). 

The issue is to show that if we have a pair {A,(j>), belonging to the class 
(|o|b,c), which solves (|ll|) on St = (0, T) x R"^ with initial data (|), and if 
the data have some additional regularity, say then this extra regularity 

persists throughout the time interval [0,T]: 

e C{[0,T],H'+'') nC^{[0,T],H''+''-^). (66) 

Now, as proved in Section 6.4], it suffices to prove this for some T > 
which depends continuously on 

Eo = Mhs + IMh.-i + WMh^ + \\Mh^-i ■ 

We shall prove this using the Picard iterates corresponding to the given data. 
It will be convenient to introduce the notation 

Ek = \\a\\Hs+k + \\b\\fjs+k-i + \\(f>Q\\fjs+k + . 

Now fix an integer K > 1, and denote by a and (3 the pointwise maxima 
of ttk and Pk, respectively, taken over all < k < K. Let us assume that the 
initial data belong to H^~^^ , that is, 

Ek < 00. 

It is proved in that for any < T < 1, there is a linear operator Wt, 
which is bounded from yj Xj [j — 1,2), and such that u ~ WtF solves 
the inhomogeneous wave equation Ou = F on (0, T) x R** with vanishing initial 
data ai t — 0. Moreover, if Ct is the maximum of the operator norms, that is, 

Ct = maxdll^Tllj;,^;,, , IIW^tII^;,^;,,) , (67) 
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then 

Ct^O as T ^ 0. (68) 

The sequence of Picard iterates {A^"^\ cj)'^"^^) is then defined inductively as 
follows. First, let A^°^ and be the solutions of DA^"^ = and D^fo) = 
with initial data (j^) , and then multiply them by a smooth bump function which 
equals 1 on the interval [0,T]. By Theorem 1], 

IIA^^-A'")!!^^ + ||AV^°^L^ < CEk, (69) 
with Ek as above. Then define 

^(m+l) ^ ^(0) ^ WTM{A^"'\(j)^"''>), 

= 0(0) i^y7V(A(™), 

Let us write 

i?!"' = II A*^^'™) II ^ + II A*^ II ^ , 

Jra) ^ |U(m) _ 11 ^ |U(m) _ , (m-1) 11 ^ 

II 1 1 ^1 1 1 ' 1 1 ^2 



Then by (|69|), (je^) and (|65|) (with = 0), we have 

4"+^' < CEo + Ct a(4'"Vo"'' > 0- 
If we choose T so small that 

2Ct a{2CEo) < 1, (70) 
then it follows by induction on m that 

< 2CEa, m>0. (71) 
Then, using the Lipschitz estimates ( p^ ) (and making a larger if necessary), 

- 2 

so the sequence of Picard iterates is Cauchy in X1XX2, and therefore converges; 
the limit is of course the unique solution {A, <j)) of our equation. 
We shall prove that, with T as in (|7C|), 



< Ck{Eo, Ek), k<K, m>0, (72) 

where Ck is some continuous function. 

Let us first see why this implies the desired conclusion ( |66| ) for k < K. The 
point is that by ([7^), the sequence of Picard iterates is bounded in the Hilbert 
space 'H^+'^fi (recall that Xi X2 = H'^'^), and therefore, some subsequence 
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converges weakly in that space. Since weak convergence in 'J-l^+'^fi implies con- 
vergence in the sense of distributions, we conclude that the strong limit {A, (f) 
agrees, as a distribution, with this weak limit. Thus, {A,(j)) belongs to 'W^^'^ , 
and this immediately gives (|66[). 

We shall prove by induction on k. 

We already have the case fc = 0, by (71). 

Now assume that k < K and that (72) holds. We claim that this implies 
(H) for fc + 1. Indeed, by (H), (|^ and (||), 

R^V" < CE,+, + CTa(4'"Vi';:l + CT/3(i?i"^). 
Taking into account ([7l|), ( [70|) and the induction hypothesis, we get 

1 p(m) l3{Ck{Eo,...,Ek)) 

R,+, <CE,+, + -R,^, + 2^(2C£;o) 

for m > 0. It now follows by induction on rn that 

using ( |69| ) for the case m — 0. 
3.2 Classical solutions 

Here we outline the proof of part (e) of Theorem ||. In view of part (d) of the 
same theorem, it suffices to prove the inductive step 

oc oc 

A, e Pi C™ ( [0, T] , =^ A,(Pe f] (7"+^ ( [0, T] , . (73) 

k=l k=l 

But since {A, cj)) solves ( |lT|) on (0, T) x R^, we have there 

a^A^ AA- 7W(yl, (/)), 

and so it is clear that (iTSf) follows from 



A,(/.G Pi c""([o,r],ff^+'=) 
fc=i 

oo 

=^ X(A,0),A/'(A,0) e p C""i([o,r],i/^+'=). (74) 



fc=i 



The key observation is of course that M and M only contain first order 
derivatives in time. Recall that M and M are sums of multilinear expressions 
in A and 4> and their first order derivatives, and terms involving Aq (</)). But 
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Aq^c/)) is determined by the elliptic equation ( ^8aD , which also contains only first 
order partial derivatives in time of (j). 

Thus, to prove (|7^), simply apply up to m — 1 time derivatives and any 
number of space derivatives, say to Ai and Af, and use the product rule 
for derivatives. It is then easy to show — we omit the details — that on each 
time-slice, the L^-norms of the resulting expressions are bounded in terms of 
(here a is a multi-index) 

ll^t^^ll^K+fc, ||^t<^||///f+'= and ||9t^9"Ao(0)||^i 

for j < 771, \a\ < K and k sufficiently large. Then one appeals to the following 
higher regularity result for Aq (0) , which is proved in section [4.3| . 

Lemma 5. Let m, M he non- negative integers. If (f> e C™+i([0, T], i/^'^), that 
is, if 

dld^<j) e C([0, T],L^) for all j <m + l and all \a\ < M + 1, 
where a is a multi-index, then 

did°AQ{(f>) eC{[0,T],H^) for all j<m and all |a| < M, 

and ||9(i9"Ao((/')||^ooQp is bounded by a continuous function of the norms 

°°([o,T]._ffA^) ^'^^ k < m + I. 

4 Elliptic estimates 

Our object here is to prove Lemmas ^ || and ||. 
4.1 Basic estimates 

We first prove existence and uniqueness for the equation 

Au-\<j>fu^~^W) (75) 



Lemma 6. Let (j) e and f & L^ . Then the equation ( [Zq ) has a unique 
(real-valued) solution u G , and 

ll^bi <2||/||i.. (76) 



Proof. Recall that , as defined in section 1.4, is a Hilbert space with inner 
product / Vu ■ Vw (by Plancherel's theorem), and that > L^. We denote 

by the corresponding real Hilbert space, with inner product / Vu • V?;. 
By definition, u £ solves (|7q) in the sense of distributions iff 



/ {Vu-\7v + \(l)\^uv)dx^ [ 'i{(j)f)vdx (77) 
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for all w G 5. Since S is dense in and 
\<j)f uv dx 
'^{(j)f)vdx < wvwj^i 



(78) 
(79) 



we conclude that u solves (|75[) iff (77) holds for all v e . Taking v — u gives 



lL2 J 



and since (0 + 6)^ < 2{a'^ + b'^) for ah a, 6 e R, we conclude that < 2A^||/||^2 
where A'' — \\u\\hi + II</'^IIl2 < 00. Therefore ( [76| ) holds, and uniqueness follows. 
Of course, u must be real, since if u solves (|7^), then solves the same equation 
with / = 0, and therefore = by what we just proved. 

To prove existence, observe that the left hand side of ( [77[ ) defines an inner 
product on 3?ij^, and in view of (|78| ), the corresponding norm is equivalent to 
the usual norm. Moreover, by (|79|), the right hand side of (^ is a bounded 
linear functional F{v) on diH^. Existence therefore follows from the Riesz rep- 
resentation theorem. □ 

Remark 3. As discussed in the introduction, our method can be modified to 
generalize the result of Cuccagna Q for MKG on M}^^ to large data in H'^ , 
s > 3/4. For this, we need the fact that ( [75| ) has a unique solution in ij^(R^) 
for (f> e H^/'^iM.^) and / e iJ-i/4(R3). A gain we multiply the equation by u 
and integrate. Using Planchcrcl's theorem we get 

\\u\\%, + \\<j)ufL2 < \\(l>u\\m/4 II/IIh-1/4 + WMlm/i ll/llff-1/4 



and since 



\\M\m/ 



< 



on R3, we get < ||0||^3/4 ||/||^ 



-1/4. 



Mm ' 

It is also easy to show that the 



operator Bq defined by (30) is bounded in for (j), Aj e H^^^ 



Next, we prove a difference estimate for ([75|). 

Lemma 7. Let (/), V' G and f,g e L^. Let u, u G be the solutions of 



Av — I ■01 V = — $J(V'5). 



Ther, 



\W-v\\Hi<U~^\\m + \\f-9\\L2 
where the suppressed constant is a polynomial in HV'lliji ^'^'^ IIsIIl^- 
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Proof. Subtracting the equations gives 

A{u -v)- 101' {u - v) = (101' - \^pf)v ~ 3 - 5)] - 3 - m ■ 
Then by a density argument as in the previous proof, 

(v(w - v) ■ V(u -v) + 101' (u - vf'j dx 

'' - \cj,f)v + 3 [0(/ - g)] + 3 [(0 - V).g]) (^x - v) dx 

< 110 - ^IIl^II'^IIl^ + II^IIlO II'^IIl^ II" - 

+ 11/ - gh. \\u - v\\^, + 110 - V^ll^. Ilsll^. \\u ~ v\\^, , 
giving the desired conclusion. □ 
We now consider the more general equation 



AM-i0ru = / 



(80) 



Lemma 8. Given G and f G L^^'^ , the equation ( |8C| ) has a unique solution 
u G , and 

Iklli/l <C||/|L4/3 (81) 

where C is independent of (j), f and u. Moreover, if 

Au-|0|'u = /, (82) 
Az;-|^|'i;-5, (83) 

where u, v,(j),ip G and f,g£ L^^^ , then 

\\u ~ v\\fj, < C(||0||^, + II^II^O ||g||^./3 110 - V'IIhi + C 11/ - 5IL4/3 , 
with the same constant C as above. 

Proof. Proceed as in the proof of Lemma ||, but with the right hand side of ( [ttI) 
replaced by — / vf dx. Thus is replaced by 



vf dx 



< V La 



r4/3 < IK'll^l 11/11^4/3 . 



Existence then follows, and any solution satisfies 



\(pu 



\l2 = - J ufdx < C||u||^i 11/11^4/3 , 



where C is independent of u, f and 0, and (^TJ) follows. 
Subtracting (||) from gives 

(A-|0|')(zi-x;) = (|0|'-|^|')z; + /-g, 
and applying (pnj gives the desired estimate. 



□ 
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Next we prove a uniqueness result in space-time: 
Lemma 9. Suppose 

(/) G C([0,r],iji) and u £ ^^([0, T], iJ^), 

and that u solves 

Au-\4>fu^0 on (0,r)xR'* 
in the sense of distributions. Then u = 0. 

Proof. Set St = (0,T) x M"* poj. every test function v{t,x) in C^{St), 

J{Vu-yv+\(j)fuv}dtdx^O. (84) 
The left hand side is a bounded linear functional in v. In fact, 



Vu • Vw dt dx 



\(f)'f' uv dt dx 



^ U\\l^(m) ll"llL?(ffi) \HLl{m) ■ 



Here we used Holder's inequality and the embedding ^ L^. 

But C^{St) is dense in L'^{[Q,T], H^), so it follows that @ must hold for 
aU i; e L^Qo, T], iji). Taking v gives 



y {|Vu|V|0nu|^}dida; = O. 



This implies Vm = 0, hence u = (i?^, as we have defined it, does not contain 
any nonzero constants). □ 

4.2 Higher regularity estimates 

Suppose 

(/.eC([0,r],iji) ^^^^ /eC([0,T],L4/3). 

By Lemma ||, the equation 

Am-|0|^u = / (85) 

has a unique solution 

ziGC([o,r],Hi). 

We shall prove the following higher regularity estimate. 
Lemma 10. Let m, M be non-negative integers. If 

dld^^ e C([0, T], iji) and dld^f € C([0, T], L^/S) 
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for all j < m and \a\ < M , then 
for j < m and \a\ < M , and 

where rja is a lower order term which depends continuously on the norms 

for all k < j and \P\ < \a\ satisfying k + \f3\ < j + \a\. Here all L°° -norms are 
taken over [0, T]. 

The proof is by induction on m and M . Denote by Pirn, M) the statement 
that the lemma holds for the pair (m, M). Since P(0, 0) is true by Lemma ^, it 
is enough, by an obvious induction, to prove 

P(m,0) => P(m+ 1,0), (86a) 
F(0,Af) =^ P(0,M+1), (86b) 
P{in+l,M),P{m,M +1) =^ P(m + 1,M+1). (86c) 

The key to proving these implications is the following: 
Lemma 11. If, for some < /i < 4, 

a^</)GC([0,T],7?i) and d^f E C{[0,T], L^/^), 
then d^u £ C{[0,T],H^) and 

II^M"llL-(ffi) ^ WfhriL^^') ll<^llLr(^^) + II^m/IIl-(l4/3) , 

where the L°° -norms are taken over [0, T]. 

Before proving Lemma [Tl], let us use it to prove (p6|). 

Proof of (|86|). Apply 9™ to both sides of (|5|). This gives 

(A-|0|^)ar«= E c,ki{dl^){dU>){d\u) + dTf. (87) 

j + k -\- I — va, I < m 

Denote the right hand side of this equation by F . By Lemma if we can show 
that F and dtF belong to C([0, T], L'*/^), then it follows that 

Sr+'w e C7([0, T],iji) 
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and 

¥T^'A\Lr(H^) < \\P\\Lri.L^r^) MLriH^) WWLriH^) + ll5*^llLr(LV3) • 

But by Holder's inequality and the embedding ^ , 

^ Irt J IIl~(li/3), 

and using the hypothesis P(m, 0) to bound we conclude that 

bounded by a continuous function of the norms ||i9t '^H^oc^^i-) 

and ||i9t /||^oc(^4/3) for all j < ™- 

Next apply dt to F. Since Z < m in (^), we can again use the hypothesis 
P{m,0) to bound ./^^ Thus, arguing as before, 

+ ll^r+ViLoo(^4/3)+i.o.t., 

where l.o.t. stands for a lower order term depending continuously on the norms 
ll^t<^IL~(iji) ^^'^ ll^t/|lL~(L4/3) for aU j < m. 

Thus P{m +1,0) holds, completing the proof of ( ^6a| ). The proof of (S6b) 
is quite similar and is omitted. 



Proof of (teed). Let |a| = M + 1. Applying 5™ 9^ to both sides of (M) gives 



(A - \cbf)drd^u = J2 cp-,5,kmd^.^mdU){d\diu) + d^d^f, m 

where the sum is over all non-negative integers j^k^l and multi-indices f3,"f,S 
such that 

l3 + '-f + 5 = a, j + k + l^m, l + \d\<m + M. 

Let F be the right hand side of (H). If F and dtF belong to C([0, T], L^/S), 
then by Lemma |ll| , 

dr+^d^u e C([0,T],iji) 

and 

||9r+'5""IL^=.(H.) ^ WnLriL^rs) Uhnm) \\dt<t>hrim) + \\dtF\\LriL^^^) ■ 



Since F and dtF are estimated just as in the proof of ( |86a| ), we will not go into 
details. The key point is that since I < m, \S\ < M+1 and /+|<^| < m+M in (p8|), 
the hypotheses P{m + 1, M) and P(to, M + 1) allow us to bound 

and||a^iaf«||^^(^,). 
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Proof of Lemma 11, Under the hypotheses of the lemma, 

d^ueL%[0,TlH'). (89) 

Before proving this, let us show that it implies the conclusion of the lemma. 

Indeed, u solves ( p5| ) in the sense of distributions on St = (0,T) x M"*, and 
if we apply 9^ to both sides, it follows that 

{A-\q^f)d^u^d^{\q^f)u + d^f. (90) 

in the sense of distributions on St- (The use of the product rule for derivatives 
is easily justified in view of (p9|).) Denote by F the right hand side of the last 
equation. Then 

and so 

^ UhriH^) II<9p'/'IIl~(hi) hhrin^) + II^m/IIl~(l4/3) • 

It then follows by Lemma ^ that the equation 

(A - \(j}\^)v ^ F (91) 
has a solution v e C([0, T], iJ^), and 

ll^llL~(iji) ^ II/IIl~(L4/3) ll'^llL~(ij"i) II^A''^llL~(ij"i) + II^A'/IIlj'°(L1/3) ■ 

Subtracting ( |9l| ) from (^ gives 

(A-|0|')(5^u-«) = 0. 

Thus d^u = w by Lemm a |9| proving the conclusion of Lemma |ri|. 

It remains to prove (pM. For technical reasons, we fix < to < 2^ and prove 
( p9[ ) with the interval [0, T] replaced by [0,to]- A similar argument works for 
the interval [to, T], giving the statement in the entire interval [0, T]. 

We shall require the following facts about the difference quotients 

j , , u{t,x + hej) — u{t,x) aO /J- _\ u{t + h, x) — u{t, x) 



Alu{t, x) = — ^ and A^u(t, x) 

where ei, . . . , 64 are the standard basis vectors of R^. 

Lemma 12. /// belongs to C([0, T], L^/'^) and the distributional derivative d^f 
belongs to C([0, T], L^/S)^ 

some < /i < 4, then 
\\Kf\\L^(l0M,L^,^)-O{l) as h^0+. 

Moreover, the same conclusion holds with f replaced by 4> and L^^^ by . 
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Proof. We have 

^lf{t) f dtfit + s) ds (L^/^-valued integral) 
h Jo 

whence 

||A^/W|L4/3< sup \\dtf{t + s)\\^,,,. 

0<s<h 

The same proof works for H^. If 1 < j < 4, then (discarding the time variable) 
Ai/(.T) = i / d,fix + se,)ds. 



Thus, by Minkowski's integral inequality and the translation invariance of the 
norm, 

||Ai/||^,/3<l|5,/||^V3. 

This is certainly valid for smooth /, and hence in general by using an approxi- 
mation of the identity. For we write 



|Ai0iii. = / icr 



e«''0 - 1 1 ^ 



h 

and note that \e'''^^ - l\ < h It follows that ||A-(0||^i < □ 
We are now ready to prove (^9|) . By the difference estimate in Lemma ^, 

l|A^u|lioo([o^to]^ij-i) ^ II/IIl==([o,t],li/3) II<^IIl-=([o,t],hi) l|A/^<?!'llioo([o^t„]^ij-i) 

for all < ft- < T — to. In view of Lemma |2|, the right hand side is 0(1) as 
h 0^ . Applying Holder's inequality in time then gives 

l|AMi.ao,t„],i^i) = as h^0+, 

so by weak compactness, there is a sequence hj —>■ such that A^ m converges 

weakly in L^([0, to], i?^) to some limit w as j — > oo. But this implies that 
AJ^ u — > V also in the sense of distributions on (0, to) x IK.*- On the other hand, 
we know that A^u — > d^u in the distributional sense as ft. — > 0, and so we 
conclude that d^u = v. This proves ( p9| ) on the interval [0,to]. 

4.3 Proofs of Lemmas |1], ^ and |5| 

First, Lemma |^ is an immediate corollary of Lemmas ^ and |^. 

Secondly, to prove Lemma ^, we apply Lemma |l^ with m = 0, AI = k and 

f = -^ 
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Since K^(f> E X2 = Ti^'^ by hypothesis, we have, in view of (p7|), 

e C([0,T],iji) and 9t5^</. e C([0, T], L^) for |a| < fc. (92) 
Thus, it suffices to check that 5^/ G C([0, T], L*/^) for |q,| < k, and 

l|5:/IL~(L4/3) < 11011^, E II^^^'^LJ' (93) 



|/3|<|a| 



where 77 is continouous. If a = 0, we have, using Holder's inequaUty and the 
embedding ^ L^, 



L2 



uniformly in < t < T, and by (17) 



\mA\dtnL2 < 



\X2 



giving ( p3| ) for a = 0. When a ^ one can apply the product rule and estimate 
each term as above. We leave the details to the interested reader. 

Finally, Lemma ^ is also proved by an application of Lemma 10, We are 
given non-negative integers m, M such that 

dld^(l)eC{[0,T],L^) for all j<m + l and aU \a\ < M + I. 

Again we set 



By Lemma |l^ it suffices to check that 

d(d^f eC{[0,T],L^^^) foraU j < m and aU |a| < M. 

Again, one simply applies the product rule for derivatives and estimates each 
term as in the proof of (p3|). 

Appendix 

Here we prove (p^. First, 

||A-i(-A)-i(H|L.(i,) < ||(-A)-i(H|L.(ij) + ||(-A)^(H|L.(^,) 



by Lemma 
both < llul 



so it suffices to show that the two terms on the right hand side are 
fjs,e [|w||^3-i,e. To this end, we apply the following theorem (stated 



here for M}'^'^ only) of Klainerman-Tataru 
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Theorem. Let 1 < p < oo, 1 < q < oo and set 7 = 2 — y — -. Assume that 

0<cr<4----, (95) 
P q 

si,S2 < 7, (96) 
si + S2 + cr = 27. (97) 

Then 



(-A)-'^/2(mw) 



where d > ^. 

Remarks. (1) When si — S2 this foUows from Theorem 4 in (see also 
Principle 3.2]). In ||^, however, the estimate was stated using the space-time 
fractional derivative operator {~At.x)~'^^^- Nevertheless, an inspection of their 
proof shows that it works equally well for (— A)"*^/^ (see Chapter 2]). In 
our statement of the theorem we have also included the end-point case due to 
Keel and Tao ||^, although we do not use this. 

(2) The asymmetric case si ^ S2 is derived as in the proof of Theorem 5 in Q. 
(The statement of that theorem contains the condition (in our notation) cr < 7, 
but an inspection of the proof shows that this is superfluous.) Let us just give 
a heuristic explanation of why the asymmetric case essentially reduces to the 
symmetric situation. Rewrite the estimate as follows: 

where D" = (— A)"/^. Denote by ^ and rj the Fourier frequencies of u and v 
corresponding to the spatial variable x. Then the frequency of the product uv 
is £^ + T], and in Fourier space, D^'^{D~''^u ■ D~''^v) is a weighted convolution, 
with weights 

1 

The idea is that the weights can be redistributed so as to get equal powers of |^| 
and \r]\. This is obviously possible if the frequencies of u and v are comparable. 
If, on the other hand, |^| » \ri\, say, then -I- ?7| ~ and so 

1 1 

< 



7 



provided si, S2 < 7 (recall that si+S2+cr = 27). Thus we are in the case a = 
and si = S2 = 7, which by Holder's inequality is reduced to a linear Strichartz 
estimate. 
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Now let p and q be defined as in section |^. Using the definition of p in (pO|), 
we see that ( p^ ) is equivalent to 1/q < {2/ 3) {6 + 2e), and the latter evidently 
holds, since 1/q < 1/4 by (|l|), ( ^9a| ) and the assumption s < 2. Thus (|4|) 
holds. 

Next we have to check that (95) holds with a — 2 (then it also holds with 



a — 3 — s, of course), but using the definition of p in 
equivalent to 

- < 26* - 1 + 4e, 

q 



we find that (95) is 



which is true by (|51|). 

Now set si = 1 + 6 and S2 = <5, where we have defined 

„ 3-cr 1 2 



(98) 



2 2p g 

With this choice, ( |97| ) clearly holds. Note that (^6|) holds provided 

1 2 

<5<1-- . 

2p q 

It suffices to check this when a = 3 — s (then it also holds for a — 2), but in 
this case it is obvious since s < 2 and, from (98), 



2p 



(99) 



It remains to check that si < s and S2 < s — 1. This is equivalent to S < s — 1, 
and again we only have to check this for cr = 3 — s, in which case it reduces to, 

by (ii, 

s 1 2 

<s-l. 

2 2p q - 



In fact. 



1 

2^ 



< s 



for by (BB), this is equivalent to 1/4 + 6'/2 + e < s/2, which holds by (49b). 
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